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The stochastic properties of fluctuating electromagnetic fields are defined in terms of 
the joint moments of the probability distribution. Their physical interpretations (coherence, 
higher order correlations) are briefly described, and the connection is indicated between 
the complete set of correlations and the quantum theory of radiation. 



1. Introduction 

This paper deals with the stochastic theory of 
radiation. It is the electromagnetic field itself 
which fluctuates here, not the medium: we are simply 
concerned with an extension to the space domain of 
the usual communication-theoretical treatment of 
signals fluctuating in time. 

As usual in a stochastic theory, we must define the 
variable and its ensemble, and describe a sequence 
of joint probability densities that specify the sta- 
tistical properties in greater and greater detail. Let 
V(x, t) represent a voltage or electric field strength 
(scalar for the time being) at point x and time t. 
Since the bandwidth of a physically realizable radia- 
tion field is never zero, V cannot be periodic in time 
but must fluctuate. Let the ensemble be a set of 
wavefields produced by one and the same source at 
different times. If ])i(Vi)dVi is the probability 
that at the space-time point (x u t{), V will have 
a value that lies between \\ and Vi+dVi, then the 
successive orders of joint probability densities are 
denned as follows: p 2 (Vi,V 2 )dVidV 2 is the joint 
probability that at the space-time point (xi, U), V 
will lie within dV\ and at (x 2 , t 2 ), within dV 2 ; 
Pz(Vi)V2,V3)dVidV 2 dV z refers analogously to three 
space-time points, etc. 

Certain weighted integrals of the joint probabilities 
often turn out to have direct physical significance: 
these are the joint moments, defined by the ensemble 
average 
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which is usually abbreviated as r U) =(F i r 2 . . . Vi). 
If the fields are stationary in time (and for conven- 
ience we shall here restrict ourselves to these), the 
products F^Ji, tx)V 2 (x2, t 2 )V-t(x 3 , t 3 ) . . . become 
V\(x u t)V 2 {x 2 , t+T)V 3 (xz, t+r') . . . , and the second- 
order joint moment, for example, is then explicitly 
written as 

r<^(F 1 (i)F 2 («+r)) 
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What is the physical significance of the several 
orders of correlation? The first order is just the 
time-average voltage (electric or magnetic) field 
at a point, which is usually zero. The second order, 
r 12 (r), traditionally called the coherence function, 
enters in the mathematical description of all inter- 
ference and diffraction effects, and of instruments 
based on these (e.g., the ordinary radio interferom- 
eter). Referred to a single point, T l2 (r) becomes 
Tn(T) = <YUt)Yi(t+ T )>, tlie autocorrelation; in 
terms of l{t), the fluctuating "intensity" (or, except 
for an admittance factor, the electromagnetic power 
flow), r u (r)=<Ii(*)>=ii, the mean (or "optical") 
intensity at the point. The third order, I heard 
recently, is being examined at present in connection 
with the correlations between incident, reflected, 
and transmitted rays in nonlinear interactions 
between light and matter. A degenerate form of the 
fourth-order correlation, applied to two rather than 
four points: 
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the cross correlation between the intensity fluctua- 
tions at two points, underlies all intensity inter - 
ferometry (Hanbury Brown-Twiss effect, etc.) and 
two-point photoelectron coincidence counting. The 
fifth order seems to have no application. The sixth, 
in a degenerate form applying to only three points, 
appears in H. Gamo's triple-correlator interferom- 
eters which, for certain specialized measurements 
(for example, the spectrum of gas discharge tubes) 
offer important advantages over ordinary amplitude 
and intensity interferometers. Orders higher than 
the sixth have not so far been found physically 
significant, but the set of all orders taken as a whole, 
i.e., the complete stochastic description of the field, 
turns out to be highly significant in the quantum 
theory of radiation. We shall return to this point 
after examining the second- and fourth-order 
correlations in somewhat more detail. 

2. Second-Order Correlations 

In figure 1, a plane, circular source of diameter 2p 
illuminates slits 1 and 2 in a screen; we will show that 
in the case of quasi-monochromatic radiation, the 
intensity I(P) at point P depends on the second-order 
correlation, the coherence function T ]2 (t). Let Yi(u) 
be a spectral amplitude component of the voltage Yi(t) 
at slit 1, and note that the corresponding power 
spectral line of Ii(t) is ii(«) = Tfi( w )*Vi( w )- The volt- 
ages at P superimpose: Y(P y u) =a l Vi(u)exp(ikr 1 ) -\- 
a 2 Y2(u>)exp(ikr 2 ) (where the a's represent the individ- 
ual slit patterns and the 1/r 2 decay), and therefore 

= \a 1 \ 2 i 1 (u) + \a 2 \ 2 i 2 (u) 

+2|a 1 ||a 2 |Re {<Fi(o>) F?(o>)>e-'-}, 

where ut is the phase difference corresponding to the 
path-length difference cr in figure 1. (All mathe- 
matical details involving the truncated functions 
have been omitted.) The total I(P) is obtained by 
integrating i(P,co) over the bandwidth Aco«co (the 
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Figure 1. Two-slit interference. 



mean frequency); noting that the last term on the 
right-hand side then represents the Fourier transform 
of the mutual power spectrum between slits 1 and 2, 
which is well known (Wiener-Khinchine theorem) to 
be the cross correlation function, we obtain 

7(P) = |a 1 | 2 / 1 + |a 2 | 2 /- 2 + 2|a 1 ||a 2 ||r 12 (0)|cosUr, 

in which the last two factors could also be written 
as Re T 12 (r). (At radio frequencies, amplitude inter- 
ferometers are used in lieu of the optical arrangement 
of figure 1: a phase shifter produces the path length 
difference, and either, as in figure 2a, a detector 
registers I(P) or else, as in figure 2b, a correlator 
consisting of multiplier and integrator reads out I\ 2 (r) 
directly.) It can be easily shown that |r 12 (0)| must 
lie between ^I 1 I 2 and 0. If we choose I 1 ^=I 2 =I J then 
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1 I(1-\-coscSt) at the upper bound and 
1 1 at the lower, as illustrated in figure 3a 
s, "complete coherence ") and figure 3d 
(no nulls, "complete incoherence"), respectively; the 
dashed curve in (a) and solid curve in (d) trace the 
individual slit pattern contained in the coefficient a. 
Cases (b), (c), and (e) are intermediate ("partial 
coherence"): for small r (near the center), the fringe 
contrast is good, but gradually it "washes out." 
The approximate path-length difference at which the 
fringes disappear for visual observation is called the 
"coherence length" cr c , where the "coherence time" 
r c ^27r/Aco; the coherence length is about 0.1 /mi for 
white light, 3 yards for a very narrow spectral line, 
200 miles for a good laser, and 2,000 miles for a well- 
stabilized klystron. 

The central fringe contrast r 12 (0), termed the 
"spatial coherence," can be measured in the Michel- 
son two-beam interferometer, figure 4a. If the source 
is incoherent (e.g., a star), the Fourier transform of 
r 12 (0) plotted against increasing separation between 
incident beams traces out the intensity distribution 
across the source (and thus also measures its diam- 



i 



i 



•**- 



I 



M 



Figure 2. Two microwave amplitude interferometers: (a) 
with phase shifter and detector D, (b) with multiplier M 
and integrator I. 
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Figure 3. Two-slit interference patterns. 
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eter). The autocorrelation rn(r), termed the "tem- 
poral coherence/' can be measured in the Michelson 
split-beam interferometer, figure 4h. Tho Fourier 
transform of V n (r) plotted against increasing path- 
length difference between minors 1 and 2 traces out 
the power spectrum of the source. 

The proof of these relations is best obtained from 
the double wave equation 

(vI-^)r 12 (r)=o, »=i,2 

which Wolf established for the coherence function; 
it is a consequence of the scalar wave equation 
satisfied by the analytic signals themselves. Thus 
the coherence function changes as light propagates 
in space: for example, complete spatial incoherence 
at the surface of a star is transformed into almost 
complete coherence by the time the light enters the 
aperture of a telescope (or else no Airy rings would 
be formed). With the wave equation one can also 
prove that the interference diagrams in figure 3 
depend on the parameter-combination Tjsp/R in 
figure 1: increasing the mean frequency of the 
incoherent source, or its diameter, or the slit separa- 
tion, or decreasing the source distance, all produce 
precisely the sequence of patterns shown in figure 3. 
(Note how the fringe contrast reappears in (e) 
after it has already been washed out in (d), though 
with less contrast than at its peak in (a), and with 
a minimum rather than maximum at the center; all 
of these effects stem from the behavior of the co- 
herence function.) 

The vector nature of electromagnetic waves is 
taken into account by writing r (2) as a matrix. For 
polarized light beams, it is sufficient to work with 
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Figure 4. Michelson interferometers: (a) two-beam, (b) split' 

beam. 



991 



where the subscripts specify which two transverse 
E or H field components are being correlated. In 
the special case [T n (0)], these four parameters 
characterize a quasi-monochromatic beam in the 
same way as the Stokes vector characterizes a 
monochromatic beam; the familiar polarization 
algebra can therefore be developed, including 
representation on a Poincare sphere, splitting an 
arbitrary beam into a fully polarized and a fully 
unpolarized part, etc. 

References to the many applications of the 
coherence calculus in diffraction, antenna patterns, 
optical imaging, radio astronomy, periodic and 
random media, etc., will be found in the last item 
listed under "Bibliography" at the end of this article. 

3. Fourth-Order Correlations 

Thermal sources consist of independently-radiating 
atoms; we know from the central limit theorem that 
all orders of the probability density distributions 
must then be Gaussian, and that all higher joint 
moments can be expressed in terms of the second. 
In particular, the degenerate two-point fourth-order 
correlation already introduced turns out to be related 
to the coherence function by 

(/i(0/2(^ + r))-/ 1 / 2 +|r 12 (r)| 2 

if the field is linearly polarized; otherwise, a factor 
multiplies the last right-hand term, for example 
0.5 if the field is unpolarized. It is clear from this 
expression that two-point intensity correlation yields 
only the magnitude of the coherence function, but 
the phase is sometimes recoverable by theoretical 
arguments, and often recoverable by additional 
measurements, for example with the three-point 
interferometer mentioned in the Introduction. 

In the Hanbury Brown and Twiss experiment, 
figure 5, |ri 2 (0)| is determined by the intensity- 
analog of the two-beam amplitude interferometer in 
figure 4a: the light or radio signals now are detected 
D and amplified A before being correlated in the 
multiplier-integrator unit. A longer baseline can 
therefore be used than alinement problems and 
atmospheric turbulence permit with the Michelson 
interferometer, but because of energy limitations 
only a handful of stars have been mapped in this 
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Figure 5. Intensity interferometer. 



way. The intensity analog of figure 4b yields 
| Tn (T) | and therefore the source spectrum (magni- 
tude only unless known to be symmetric) . 

Of the several other experiments that depend on 
intensity correlation, I will mention only that of 
Alford and Gold, in which power spectrum modula- 
tions are observed when a beam is recombined after 
a path-length difference in excess of the coherence 
length. 

The fourth-order correlation of nonthermal sources, 
such as lasers and klystrons, does not depend on the 
coherence function; in fact, an ideal single-mode am- 
plitude-stabilized oscillator should produce intensity 
fluctuations of zero correlation. 

In optics, two-point coincidence counting is often 
used in lieu of intensity correlation (Pound and 
Rebka experiment, etc.). The equation that estab- 
lishes a stochastic relationship between the proba- 
bility distribution of the photoelectrons in time and 
the intensity fluctuations is 

p(t)dt=al(t)dt, 

where a depends on the efficiency of the photode- 
tector. Originally suggested by Purcell on a semi- 
classical basis applicable to monochromatic fields 
only, this relation is now known to hold for poly- 
chromatic fields within the accuracy of first- 
order quantum-mechanical perturbation theory. A 
straightforward calculation shows that the prob- 
ability of finding n photons in a time T consists, 
for thermal radiation, of two terms: a classical par- 
ticlelike Poisson distribution added to a wavelike 
Bose-Einstein distribution (not to be confused with 
the basic Bose-Einstein statistics which underlies all 
photon distributions). For ideal laser or klystron 
radiation, the distribution is completely Poisson, 
i.e., like shot noise — which explains why there can 
be no two-point intensity correlation in this case. 

Returning to thermal radiation, the variance in 
the number n T of photoelectrons ejected during T 
follows from the distribution 



(An T ) 2 =n T +a 2 (I T T) 2 , 

which again is the sum of a term due to classical 
particles and one due to classical waves — just like 
Einstein's celebrated blackbody fluctuation formula, 
but now valid also for radiation that is not in thermal 
equilibrium. It can be shown from this expression 
that the two-point coincidence counts, for T^>^>t c 
and linear polarization, are given by 

An 1 An 2 =a 2 r c T\T 12 (0)\ 2 

(with a factor 0.5 for unpolarized light). We have 
thus available a third method for measuring the 
intensity distribution across incoherent sources; a 
related procedure yields spectral information. 

The fourth-order correlation propagates through 
space in accordance with a four-fold wave equation. 
(And it is true in general that an wth-order correla- 
tion satisfies an n-fold set of wave equations.) 
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4. Connection With the Quantum 
Theory of Radiation 

It is well known that statistical mechanics, which 
treats ensembles of classical particles, contains more 
physics than one might expect from an n-body prob- 
lem in classical mechanics; in particular, it encom- 
passes irreversibility, whereas the equations of 
classical mechanics are reversible. This enrichment 
in physical content comes from the assignment of 
probability distributions to quantities that are 
sharply defined classically (for example, the velocity 
of molecules), and from certain auxiliary concepts, 
especially equip artition and the rules for counting 
degrees of freedom. 

Here, too, we have "enriched" a classical theory 
with probability distributions; can auxiliary concepts 
now be added so as to produce a physical theory of 
larger scope than classical electrodynamics? The 
answer is affirmative: the stochastic theory of classi- 
cal radiation we have just sketched can be further 
developed so as to encompass quantum effects; in 
fact, recent work by Glauber and Sudarshan shows 
that, for all linear interactions, the enlarged sto- 
chastic theory is isomorphic with the quantum theory 
of electromagnetic fields. 

The basic correspondence is found to be that be- 
tween the complete sequence p iy p 2} pz • • ■ of prob- 
ability distributions with which we introduced the 
classical ensemble, and the quantum- mechanical 
density matrix; the classical correlations are then the 
expectation values of the quantum-mechanical ob- 
servables. Two auxiliary concepts are implied by 
this equivalence: First, because the quantum-me- 



chanical operator representing the electromagnetic 
field is complex, use of the analytic signal becomes a 
necessity rather than a convenience. And second, 
the Hermiticity of the density matrix forces the 
probability distribution to become negative over cer- 
tain ranges of the field variables; these turn out to 
be unobservable because they violate the uncertainty 
principle. 

Could it be that the confrontation of the engineer 
with areas of knowledge previously reserved to the 
pure physicist, which occurred during the early 
'Fifties in solid-state electronics, and more recently 
in maser and laser technology, will now repeat itself 
in the field of radio and optical wave propagation? 
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